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?—i ' Abstract 



We propose a consistent setup for a holographic dual of Bjorken flow of strongly 
coupled large- N c Af=4 SYM-theory plasma. We employ Eddington-Finkelstein 
type coordinates for the dual geometry, and we propose a late-time expansion 
there. We construct the dual geometry order by order, and we show that the 
transport coefficients are determined by the regularity of the geometry. We also 
show that the dual geometry has an apparent horizon hence an event horizon, 
which covers the singularity at the origin. We prove that the dual geometry is 
regular at all orders under an appropriate choice of the transport coefficients. 
Our model is a concrete well-defined example of time-dependent AdS/CFT. 
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1 Introduction 

The application of the Anti-de-Sitter space/ Conformal Field Theory (AdS/CFT) corre- 
spondence [U [2] to physics of quark-gluon plasma (QGP) has been one of the active re- 
search fields after the observation of strongly interacting QGP at the Relativistic Heavy- Ion 
Collider (RHIC). Especially, holographic description of time-dependent systems is impor- 
tant, because the QGP at RHIC and the one expected to be observed at the Large Hadron 
Collider (LHC) are time-dependent. If we attempt to simulate the time-dependent QGP 
by using lattice QCD, we need a huge computational power. A useful framework for the 
time-dependent QGP is a relativistic hydrodynamics which describes the macroscopic na- 
ture of the plasma. However, microscopic information on the plasma, such as equation 
of state and transport coefficients, has to be provided by separate computations based on 
the microscopic theory. The AdS/CFT correspondence is an interesting research field for 
the time-dependent plasma, since it describes both the microscopic and the macroscopic 
nature of YM theories within a single framework. 

An important series of works along this direction is studies on a Yang-Mills (YM) theory 
fluid which undergoes the Bjorken's boost-invariant one- dimensional expansion (Bjorken 
flow) [3], initiated by Janik and Peschanski [4j. The Bjorken flow is a standard, simplest 
model which well captures the nature of QGP. Construction of a time-dependent framework 
of AdS/CFT is a challenging task in itself, and the Bjorken flow of an M = 4 SYM 
fluid has been studied with a late-time (large proper-time) approximation. A way of 
taking the late-time limit in the dual geometry has been proposed in Ref. [4J and a proper 
late-time expansion has been found to be given with respect to r -2 / 3 , with r being the 
proper-time [5]. The analyses have been done to the zeroth ordeijj in the expansion [H 
E], to the first order [5], to the second order [7J, and to the third order [SJ E] - It has 
been proposed in Ref. [4] that the regularity of the holographic geometry determines the 
hydrodynamic parameters of the corresponding gauge-theory fluid. Indeed, the equation of 
state [I], the shear viscosity [7J and the relaxation time [8] have been uniquely determined 
from the regularity. Amazingly, the results agree with those computed by other methods. 
(See Ref. [TDj and reviews [TT] for the shear viscosity, and Refs. [TS1 [T31 dH [T5] for the 
relaxation timaj.) The expectation value of the field strength squared (TrF 2 ) has also 
been determined to be zero at the leading order [6| and to the second order (more precisely, 
to the order of r~ 3 ) |9] by the regularity. See Refs. [17] for other related works. 

Despite the above success, there are still issues what we need to understand better: 



1 We call the leading order "zeroth order." Then, our first-order (second-order) geometry describes the 
first-order (second-order) dissipative hydrodynamics. 

2 A new second-order hydrodynamics has been proposed in [TU Q~5] , while the hydrodynamics in Ref. [H] 
is based on the Miiller-Israel-Stewart theory [16] (see also Ref. [13] for some comment). However, as far as 
the obtained values of the transport coefficients are concerned, the difference among them merely comes 
from the definition of the coefficients, and their results are consistent with each other. For example, one 
of the new transport coefficients, which is denoted by Ai in Ref. [14] . has been determined by combining 
the result of Ref. [5] and an independent computation of the relaxation time [TJ] . The result agrees with 
that of Ref. [15]. 



1) In the previous works [U [6j El 13 EJ [9] (see also related works [17]). the physical interpre- 
tation has been made by assuming that the background geometries are time-dependent 
black holes. However, it is non-trivial to see the presence of an event horizon in a time- 
dependent geometry. As far as the authors know, any rigorous proof of the presence 
of the event horizon on the dual time-dependent geometry has not yet been reportedo 
We need to examine whether they are really black holes or not very carefully. 

2) One natural question is why the regularity of the holographic geometries determines the 
physical values of the parameters. From the string-theory point of view, the presence 
of the singularity merely means a break-down of the supergravity approximation, and 
does not necessarily mean the break-down of the physics itself. 

3) It has been claimed in Ref. [9] that there is a logarithmic singularity at the third order 
in the dual geometry which cannot be removed within the framework of 10-dimensional 
I IB supergravity. 

In the present paper, we solve the problems 1) and 3), and discuss a possible idea which 
may answer the question 2). A key object in the present paper is an apparent horizon, 
that is defined as a boundary between trapped and un-trapped regions. For the problem 
1), what we need to show is the presence of the event horizon in the given background. 
However, analysis of event horizon is hard in a time-dependent geometry in general since 
event horizon is defined globally. A more convenient object is the apparent horizon which 
is defined locally. In this paper, we compute the location of the apparent horizon explicitly 
in a newly proposed dual geometry, and we show its presence. Since the presence of an 
apparent horizon is a sufficient condition for the presence of an event horizon [20], we prove 
that the dual geometry is really a dynamical black hole. 

For the problem 2), we point out that the idea of cosmic censorship hypothesis (CCH) 
may be helpful. The cosmic censorship hypothesis [211 [22] says that naked singularities do 
not appear in any physical process in the gravity theory, and all the singularities which are 
created in the dynamical process must to be "hidden" by the event horizon. Although CCH 
is a conjecture, no definite counter-examples of CCH in asymptotically AdS spacetimes has 
been found so far. (See, for example, Refs. [21] ■) The reason why the regularity of the 
dual geometry can be a physical condition becomes clear if CCH holds. Suppose that we 
choose a certain value of a hydrodynamic parameter and we find the corresponding dual 
geometry has a naked singularity. CCH says that such a geometry cannot be created by 
any physical process in the gravity side. This means that we have no way to create such a 
plasma with that particular value of the parameter as a result of any physical process of 
the YM theory as far as the duality holds. This explains why that value of the parameter 



3 A "local event horizon" has been defined in Ref. [TS] and studied in Refs. [181 IT9] . 

4 For example, we have a naked singularity in the Constable-Myers (CM) geometry [23] which is static. 
Here we interpret the statement of CCH that "any naked singularity is not created dynamically by starting 
with a regular geometry," and we do not take the CM geometry as a counter-example of what the CCH 
means in the present paper. 



has to be discarded. Then, precise examination of the location of the event horizon is very 
important to judge whether the singularity is covered by the horizon or not. Since the 
location of the event horizon in a time-dependent geometry is non-trivial, we need careful 
examinations. Again, the apparent horizon gives important information. The location 
of the apparent horizon provides a bound for the position of the event horizon since the 
apparent horizon exists necessarily inside (or on top of) the event horizon [20] . 

In the present work, we find that the Fefferman-Graham coordinates which have been 
exclusively utilized for the holographic dual of Bjorken flow are not appropriate for the 
description of the apparent horizon. Then we propose to construct a dual geometry on 
ingoing Eddington-Finkelstein type coordinatesjj We define a late-time approximation and 
we explicitly construct the dual geometry to the second order of the late-time expansion. 
We find that the regularity of the bulk geometry except at the origin uniquely determines 
the transport coefficients. One technically new point is that the transport coefficients are 
determined from the regularity at one order lower than those in the previous works [7J [8] . 
A key quantity for doing this is a Riemann tensor projected onto a regular orthonormal 
basis. This brings a technical benefit for computations of transport coefficients of more 
complicated models whose higher-order geometries are not easily obtainablejj Furthermore, 
we prove that the dual geometry is regular (except at the origin) for all orders if we choose 
the transport coefficients appropriately. We show that such a choice exists at every order. 
Therefore the logarithmic singularity pointed out in Ref. [9] is absent from the newly 
proposed geometry, hence the problem 3) is solved. Our interpretation is that the late- 
time expansion on the Fefferman-Graham coordinates is ill-defined. We also show the 
presence of the apparent horizon hence the event horizon which covers the singularity at 
the origin. We compute the location of the apparent horizon explicitly to the second order. 
We show that our geometry is non-static even from the viewpoint of local geometry. We 
conclude that our geometry is a dynamical black hole and the present model is a concrete 
well-defined example of time-dependent AdS/CFT. 

Before closing the introduction, let us specify the assumptions we shall use in this 
paper. Let (r, y, x 2 , x 3 ) be the local rest frame (LRF) (the comoving frame) of our Bjorken 
flow on which the fluid is at rest. Here, r is the proper-time, y is the rapidity, x 2 and 
x 3 are the perpendicular directions to the collisional axis. (See also Appendix O) We 
assume that the fluid extends homogeneously in the perpendicular directions, and we have 
translational and the rotational symmetries on the (x 2 ,x 3 )-plane. Another assumption is 
the presence of the boost invariance which is the translational symmetry in the y direction. 
In the realistic QGP, the boost invariance is realized at the central rapidity region where 
y is small. However, we assume that the boost invariance holds in the entire region of y 
in our setup for simplicity. The symmetries on the (x 2 , x 3 )-plane are also approximately 
realized at the vicinity of the collisional axis in the central rapidity region (namely, at the 



5 Ingoing Eddington-Finkelstein coordinates are used in Refs. [15l [T8l [T9] to construct a holographic 
dual of a plasma perturbed around a static configuration. See also Refs. [25JI26]. 

6 For example, the shear viscosity of the Bjorken fluid at finite coupling may be obtainable more easily 
than the work of Ref. [27]. (See also Ref. [28].) 



central part of the plasma) in the case of the central collision. It should be understood that 
we are investigating the nature of the fluid in this region, if one attempts to compare with 
the realistic QGP. We also assume that the expansion rate of the fluid is slow enough so 
that the hydrodynamic description is valid. In other words, we assume the presence of the 
local thermal equilibrium. Of course, our system is time dependent and dissipative; our 
system is not at the thermal equilibrium, although all the portions of the fluid share the 
same (time-dependent) temperature because of the symmetries we have assumed. Here, 
the "local thermal equilibrium" means that the expansion rate of the fluid is slow enough 
comparing to the typical microscopic time scale of the fluid (say, the relaxation time). 
Since the expansion rate of the Bjorken flow becomes slower and slower along the time 
evolution, we assume that r is large enough comparing to the microscopic time scale. 

The organization of the present paper is as follows. In Section [2J we point out the 
difficulties in the Fefferman- Graham coordinates. In Section [31 we summarize how the 
hydrodynamic equation and the equation of state are encoded in the Einstein's equation 
and the bulk theory. In Section HJ we propose a new recipe to construct a dual geometry 
on Eddington-Finkelstein type coordinates based on the late-time approximation. If we 
parametrize the dual geometry naively, it is not manifest how the dual of empty fluid is 
reduced to pure AdS. We propose a parametrization which makes the reduction manifest. 
In Section [5j we construct the dual geometry explicitly and analyze to the second order. 
The regularity of the geometry for all orders is discussed in Section E] In Section [7J, we 
compute the location of the apparent horizon and prove the presence of the event horizon. 
The non-staticity of our geometry is briefly commented in Section [HJ We conclude in 
the last section. A number of overviews that may be useful for the readers are given in 
Appendix. 



2 Problems in Fefferman- Graham coordinates 

The Fefferman-Graham (FG) coordinates are very useful for the description of the holo- 
graphic renormalization [29J , and have been used to describe a holographic dual of Bjorken 
flow [H El El [7J, [8J . (See also Refs. [T7J.) However, we point out that there is a crucial prob- 
lem in FG coordinates which prevents us from investigating (dynamical) apparent horizons 
in the dual geometries: we cannot see any trapped region on this coordinates. 

The apparent horizon is defined as the boundary between the trapped and un-trapped 
regions. An intuitive but not very rigorous explanation is as follows. The trapped region 
is the region where the light emitted outwards propagates inwards due to the gravitational 
effect of the black hole, while the un-trapped region is where the light emitted outwards 
propagates outwards. The apparent horizon is the boundary between the two regions. 

To examine the location of the apparent horizon, it is convenient to consider null normal 
expansions 9± and their product B = e^(9 + 6L, where e^ is an appropriate normalization. 
The expansions 6± are defined by 9± = C± log ft where /i is the unit volume of the spatial 



3-surface and C± are the Lie derivatives along the null normal vectors. (See Section 17.11 
for the details and the precise definitions.) An intuitive but not very precise explanation 
is as follows. Let us consider a surface which is "perpendicular" to the light array. 9 
measures how the volume of the surface grows along the propagation of the light. The 
light propagates outwards if 9 > 0, while it propagates inwards if 9 < 0. Then, the trapped 
region is defined as the region of > where the light propagates inwardsj regardless of 
their emitted directions (±). The un-trapped region is defined as the region of < where 
the light emitted outwards propagates outwards and the one emitted inwards propagates 
inwards (as they do on the flat spacetime). The location of the apparent horizon is given 
by = 0. 

To demonstrate the problem of FG coordinates, let us consider the geometry given 
in Ref. [I] and examine whether it has an apparent horizon or not. What the authors of 
Ref. [1] have found is the following. Suppose that the proper-time dependence of the energy 
density of the Bjorken fluid werqj e(r) ~ r~ l . The dual geometry in the large proper-time 
region is obtained to be [I] 



ds 2 



gl^(r,z)dx fl dx u + dz 2 



with g™ 



diag(-e<\ tV, e c ~, e c ) 



(1) 



where z is the 5th coordinate and 



— log 
25 S 



5v 4 ) 



4U+5 



:i + 5v*y- 5 



' b= 25 l °' 



Sv 



4\l-l+8 



[l + Sv 



4\l-l-5 



— log 
25 S 



1 _ Sv ^-l+l/2+S 

(1 + 5v 4 )- 1+l / 2 - s 
(2) 



Here v = zt~ 1 ^ 4 and 5 = a/(3/ 2 — 8/ + 8)/24. They found that the regularity of the 
Riemann-tensor squared (which we call Kretschmann scalar in this paper) singles out the 
correct physical value I = 4/3. 

Let us attempt to compute the location of the apparent horizon. for this geometry 
is given byj 



= - 



9 



2„,8 



(^ 4 -3) 2 + |(/-4/3) 2 t; 
9 - v s - |(/ - 4/3)V* 



+ 0{r—). 



(3) 



Interestingly, at the leading order can be zero only when I = 4/3. If we set / = 4/3, 



= - 



9 



(3 ~ v 4 ) 

3 + v 4 



-i 2 



+ 0(r—). 



(4) 



7 We do not consider white holes here. 

8 Of course, we know that I has to be 4/3 as a consequence of the equation of state and the hydrodynamic 
equations of the perfect fluid [3]. (See for example, the first term of (|114|l .) 
9 Here, we assume that < I < 4 because of the positive-energy condition [4]. 



and the candidate for the position of the apparent horizon is v = z/t 1 ^ 3 = 3 1 / 4 ; one 
may conclude that the presence of apparent horizon singles out the correct proper-time 
dependence of the energy density (hence the correct equation of state). However, we cannot 
conclude at this stage. One should notice that the G in (jl]) is always negative or zero: there 
is no trapped region. 

The origin of the problem we have encountered above is understood by considering both 
the static AdS black hole (AdS-BH) on FG coordinates and that on the Schwarzschild-type 
coordinates. A metric of a static AdS-BH on FG coordinates is given by 

i f (i-(zM,) 4 ) 2 

* 2 1 1 + (z/zq) 



ds = -5 1 - -, , , \1 dt + (! + {*/**) W +dz 2 }, (5) 



where Zq is the location of the event horizon. We can switch to the Schwarzschild-type 
coordinates through the coordinate transformation, 

r- 1 = - = (6) 

VI + (*/*b) 4 

The resultant metric is 

d* 2 = ~r 2 [1 - (ro/r) 4 ] dt* + rW + ^ fl _ ( ' o/r)4] ^ (7) 

The important point is that 

rz = ^(z/zo) 2 + (zq/z) 2 > V2, (8) 

and the equality holds at the event horizon: the entire region of the ^-coordinate covers 
only outside the event horizon (namely, only the un-trapped region) in the Schwarzschild- 
type coordinates. The points at z and Zq/z on the FG coordinates are mapped to the 
same point on the Schwarzschild-type coordinates outside the horizon. In the dynamical 
setups, the map between FG coordinates and the Schwarzschild-type coordinates are more 
complicated. However, we have seen explicitly the same problem in the dynamical example 
above. 

Let us go back to the geometry (J2J). Since we cannot show the presence of the trapped 
region, the point which satisfies = is only an candidate for the location of the apparent 
horizon; we need to postpone the conclusion until we show the presence of the trapped 
region. Furthermore, we cannot conclude the absence of the apparent horizon at / ^ 4/3, 
since we have not examined the entire region of the full geometry. 

One may expect that the trapped region may appear if we include the higher-order 
contributions of the late-time expansion. However, we find that the late-time expansion 
fails at the vicinity of v — 3 1 / 4 on this foliation and we cannot examine the location of the 
horizon in a well-defined way. See for the details, Appendix |A] 



3 Construction of dual geometry 

In the previous section (and in Appendix[A}, we have observed difficulties in the holographic 
dual of Bjorken flow on the FG coordinates. Then, we need to construct the dual geometry 
on a better coordinate system based on a well-defined approximation. We propose to 
construct the dual geometry on the Eddington-Finkelstein type coordinates, where the 
trapped region and the untrapped region are packed into a single coordinate patch. In this 
section, we summarize how to construct the dual geometry from the boundary data. 

The dual geometry has to be a solution to the lOd type IIB super-gravity equation. 
However, for the systems we consider^! the super-gravity equation is reduced to a five 
dimensional (5d) Einstein's equation with a negative cosmological constant A = —6 [29] n l: 

K - \fiR - %, M = o. (9) 

To fix the geometry, we need to specify the boundary condition which is the input for our 
theory. Let us clarify our working standpoint about what our inputs are. We take the 
Bjorken flow as an input of the theory in the present paper. The flow is specified by using 
the local rest frame (LRF) (the comoving frame) on which the fluid is at rest. The LRF 
metric for the Bjorken flow is given byt 12 l 

dsj = -dr 2 + r 2 dy 2 + dx 2 ± , (10) 

where r is the proper-time, y is the rapidity and x± denotes the perpendicular directions to 
the collisional axis. The LRF, (r, y, x 2 , x 3 ) in our convention, is given by the boost transfor- 
mation from the cartesian coordinates, and their relationship is {t,x l ) = (r coshy, r sinhy). 
(See for the details, Appendix O) Since the Bjorken flow has translational invariance in 
y direction (which is called boost invariance), we assume that our metric does not depend 
of y. We also assume that the fluid is homogeneous in x± directions, hence the metric 
is independent on x±, too: we assume that the metric is a function of only the time-like 
coordinate and the radial coordinate in the bulk. 

Based on the above standpoint, the boundary condition is given by the LRF (fTOj) ; we 
take the Dirichlet boundary condition for the metric rather than the Neumann bound- 
ary condition. In this case, it is known that we need to add the Gibbons-Hawking- York 
boundary term to the bulk Einstein-Hilbert action to make the variational principle well- 
defined EH. 



The precise dictionary between the boundary condition and the 4d geometry is given 
by the GKP-Witten prescription [2J, where the non-normalizable mode of the bulk metric 



10 We assume that the dilaton and the RR 5-form field strength do not depend on time and they are 
the same as the static case. They solve the super-gravity equation as far as ((9]) is satisfied. 

11 Here, A = — d(d — 1)/(2Z§) , where d = 4 is the dimension of the boundary theory and Iq is the length 
scale of the geometry. We set Iq = 1 in the present paper. We take the convention of the curvature tensor 
in such a way that R < for AdS. 

12 This is the Rindler spacetime in general relativity language. 



is identified with the 4d metric. To be specific, let us consider an Eddington-Finkelstein 
type coordinate system: 

ds 2 = r 2 g T+T+ drl + 2dr + dr + r 2 g yy dy 2 + r 2 g x±x± dx]_, (11) 

where g^ are functions of only the time-like coordinate r + and the radial coordinate r. 
Then, 

9ij\r—>oc v / 



is identified with our 4d LRF metric (TTOl) . Once the boundary metric is given, the expec- 
tation value of its conjugate quantity, the 4d stress tensor, is obtained by differentiating 
the bulk action with respect to the boundary metric. We shall see in Sections |5] and [6] that 
the stress tensor is indeed determined uniquely (up to overall normalization). 

The differentiation of the bulk action with respect to the boundary metric is considered 
in a covariant way in Ref. [30]. Let us introduce a regularized boundary which is a constant- 
r surface, and we define the induced metric on the regularized boundary as %, ,,, (which 
contains r dependence). The covariant dictionary obtained in Ref. ^0\ is then 13 ! 



1 » V \ 4^2 ' ; 



K^y - K^^y - 37 Mi , + -G^y 



(13) 



where G^ v is the boundary Einstein tensor (with zero cosmological constant) with respect 
to 7 M j,. Kfj, u is the boundary extrinsic curvature which is defined as 



1 

2 ( 



Kp, = ~-(V M n„ + V,fv), (14) 



where V M is the covariant derivative with respect to 7^ and n^ is the outward-pointing 
unit normal vector to the regularized boundary. See Appendix [B] for more details. The 
first two terms in the bracket of (TT3|) came from the Gibbons-Hawking- York boundary term 
[3T] while the last two terms are the counter terms which have been introduced to remove 
the divergence [30]. The finite contribution from the counter terms is crucial to get the 
correct result. Eq. (TTBl is essentially given by the normalizable mode of the bulk metric. 
The precise map between the asymptotic behaviour of the bulk metric and the stress tensor 
shall be given in Section 13.21 

The remaining task is to interpret the resultant stress tensor in terms of the hydrody- 
namics. The interpretation is given by comparing the stress tensor with its hydrodynamic 
definition. From the hydrodynamic computations, we obtain the following result for the 



13 We have inserted r 2 in order to define the 4d stress tensor in an r-independent way. We have also 
used the relation between the 5d Newton's constant G5 and N c , which is given by (87rG5) _1 = N 2 /(Att 2 ) 
in our convention. See also Refs. [29), [32] for the dictionary on the FG coordinates. 



Bjorken flow: 

T TT /e = r- 4 / 3 J2 4 ] r- 2k/3 = r^ 3 - l^ + egV 8 ' 8 + • • • , (15) 



fc>0 



T yy /e Q = '-r^ - 2 Vo + ^r^ + ■ ■ ■ , (16) 

T x±x J^ = \t^-\^t-^ + ..., (17) 



where we identify r + = r at the boundary. Here, 



e (2) 

6 



H 2 + 4A 
, A = A 1 — ri T n , 



and €q is the overall normalization of the energy density of the fluid. r] , r^, A 5 are the pa- 
rameters which are proportional to the shear viscosity, the relaxation time and a transport 
coefficient introduced in Refs. [THUS], respectively See, for more details, AppendixO The 
comparison between ([13]) and ([15]) . (fT6l) . ffTTl) enables us to read the transport coefficients. 

It is important to realize that we have not introduced equation of state nor hydrody- 
namic equation into the gravity-dual side by hand. They are automatically encoded in the 
dual theory. We shall demonstrate this in the next subsections. 



3.1 Hydrodynamic equation from Einstein's equation 



As is advertised in Ref. [To], the hydrodynamic equation is given as a consequence of 
the Einstein's equation in the gravity dualo We present a general derivation of the 
hydrodynamic equation here. 

We point out that the 4d stress tensor is related to the 5d Weyl tensor] 1 ^ The Gauss 
equation for the r-constant surface gives 

Ra^u = (5) R K X P ela K JpW + K^Kfr ~ K av K^. (19) 

By contracting this equation and by using the bulk Einstein's equation ^R^ — h^Rg^ — 

[1U 



6gn U = 0, we obtain the following relationship [3 /1 '"' ! ' 



G> = 3 lfiu + K^K - K^K v a + i (K al3 K a P - K 2 ) liiV - (5) C^ ul3 h a h^ (20) 



14 See Refs. [33] where the dynamics of the fluid is obtained as a consequence of the Einstein's equation. 

15 In this section, we put ( 5 ) for the 5d geometrical quantities (defined with respect to g^ v ) to avoid 
confusion. The quantities without ( 5 ) should be understood as the 4d quantities which are defined with 
respect to 7^, in this section. The quantities without ' 5 ' in other sections are five-dimensional ones, for 
notational simplicity, if it is not specified. 

16 See also Ref. (35;. 
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where ^'C IMM/ is the 5d Weyl tensor. Let us define K^ = K fJjU +'j fll/ for convenience. Then 
the above equation can be rewritten as 

G^ + 2K^ V - 2Kj^ - 67^ 

= k^k - k m k v a + 1 (k a pk af3 - k 2 } liiv - ^c m ^h a h^. (21) 

Now the boundary condition at r — > oo yields K M! , = 0. Hence we have the final expression 



■^■/xv **-iiiu ^ifiu T" r) '- T /ii/ 



1 
-I 

2 



■i (5) C MQ ^n Q ^. (22) 



The left-hand side is nothing but the 4d stress tensor and it is now given by the projected 
5d Weyl tensor. It is obvious that the right-hand side satisfies the traceless condition 
(which gives the equation of state in terms of hydrodynamics) because of the traceless 
property of the Weyl tensor. 

Furthermore, the Codazzi equation is 

V,K\ - V U K = ^R a ^ a u n^. (23) 

By using the bulk Einstein's equation again, we easily find V pK^v — V^if = for any 
r-constant surface. The Bianchi identity for the 4d Einstein tensor G^ v and the Codazzi 
equation yield the conservation law (the hydrodynamic equation) for (122|) . Thus, we have 
shown that the equation of state and the hydrodynamic equation are given as a consequence 
of the 5d Einstein's equation. 

3.2 Stress tensor from asymptotic geometry 

We demonstrate, based on a concrete example, that the hydrodynamic equation and the 
equation of state are obtained by solving the Einstein's equation at the vicinity of the 
boundary. Let us expand g T+T+ with respect to 1/r: 

- g T+T+ ( r+ , r) = 1 + a^r- 1 + ■■■ + a (4 V" 4 + • • • . (24) 

We also expand g yy and g x±x± with respect to 1/r and substitute them to the Einstein's 
equation. The Einstein's equation relates the coefficients of the expansions, and we found 
that they are written as 

(r +1 r) = l + a^r- x +(^P^-d r+ a^)r- 2 + a^r- 4 + 0(r- 5 ), 



9t + t + \ 

9yy( T +,r) 



4 
r 2 + T + (r + aW + 2)r~ l + - A (r + a^ + 2) 2 r" 2 



l««l'_, l/ (4) + 3 ( 



9* ± .Ar+,r) = l + a^r- 1 + K --^r- 2 --[a^ + ^d T+ a^)r + )r- 4 +0(r- 5 )(25) 
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Notice that a^ and a^ depend on r+. 

Substituting the above expressions into (Tl3l) . we find 



T 

J. TT 


= 


-Ec 


» w , 


T 
± yy 


— 


Et 2 


(4) + Tdra W 


T 


— 


-E 






E 


3N* 



(26) 

where 

3 /v 2 

2 4,- < 27 » 

We have rewrote r + as r, since r + at the boundary is identified with the proper-time of 
the fluid. Notice that aP~' does not appear in the stress tensor. We shall see that a^ 1 ' 
corresponds to a gauge degree of freedom, in Section 

It is interesting that all the components of the stress tensor are given by using only 
a^ . The relationship among the above three components of the stress tensor agrees with 
the one given in Eq. (5) in Ref. [3]. We can show in general that T yy and T x±x± of the 
Bjorken flow are expressed by using T TT as f)26p . if the stress tensor is conserved and 
traceless. Therefore, the above is a concrete manifestation of what we have concluded in 
the previous subsection: the Einstein's equation at the vicinity of the boundary yields the 
hydrodynamic equation together with the equation of state. 

At this stage a^ is a function of r which cannot be determined from the LRF and the 
symmetry of the system. However, it is determined from the regularity of the geometry as 
we shall see in Section [5] and in Section [6j In order to identify a^{r) with (TT5l) . we need 
to define a well-defined 1/t + expansion (which we call the late-time expansion) in the 
bulk theory. We shall propose the late-time expansion on the Eddington-Finkelstein type 
coordinates in the next section. 

4 Our proposal: Gravity dual of Bjorken flow on 
Eddington-Finkelstein type coordinates 

We propose a late-time expansion on the Eddington-Finkelstein type coordinates in this 
section. We also summarize our proposal on the construction of the dual geometry. 

4.1 Basic philosophy 

A good starting point for us is to re-interpret the work of Ref. [1] in the following way. 
The static AdS-BH on FG coordinates, given by (|3J), has a Hawking temperature T# = 

12 



v2/(7TZo). On the other hand, we know that the temperature of the Bjorken fluid depends 
on the proper-time as T ~ r -1 / 3 |3J. Then, the dual geometry of the Bjorken fluid may 
be described by replacing z with u^V 1 / 3 at (JSj), where w is a constant. However, this is 
not enough since the boundary coordinates should be the LRF. The Minkowski metric on 
the local rest frame is given by (1101) . Therefore the dual geometry may be given by 



ds 2 



1 " <Sr 



w, 



T 4/3 



-dr 2 +[l + w * 



4/3 



T 



{r 2 dy 2 + dx 2 ± ) + dz 2 



m 



Indeed this is what the authors of Ref. [3] have obtained. The energy density of the fluid 
is proportional to Wq. We can easily see that ([28]) is reduced to a pure AdS geometry if we 
take limit of Wq — > 0. This is consistent with the picture that the fluid becomes empty at 
this limit. 

Let us follow the same procedure on the Eddington-Finkelstein coordinates. A static 
AdS-BH on the ingoing Eddington-Finkelstein coordinates is given by 



ds 2 



r ) 



dt 2 , + 2dt+dr + r 2 dx 2 



(29) 



where t + is the time-like coordinate and r is the radial coordinate. The Hawking temper- 

-1 to 

ature is given by Th = tq/tt. Let us replace r with wr + where w is a constant and we 
regard the boundary value of r + as the proper-time. We also replace the boundary metric 
with that of the LRF. We reach 



ds 2 = -r 2 



1- 



w 



rr 



1/3 



dr\ + 2dr + dr + r 2 {T 2 ,dy 2 + dx] 



(30) 



as a candidate for the dual geometry. A natural interpretation is that w A is proportional 
to the energy density of the fluid. 



One may notice that (130]) does not reach pure AdS geometry at the limit of w = 0. 
A crucial difference between (130]) and (128]) is the presence of the off-diagonal component 
2dr + dr in (130]) which mixes the time-like coordinate and the radial coordinate. We may 
improve (I3"U]) by modifying the (y, y) component: 



ds 2 



w 



rr. 



1/3 



dr\ + 2dr + dr + r 2 r\ (l-\ j dy 2 + r 2 dx] 



(31) 



Then (I3"T]) is reduced to an exact pure AdS geometry at w — > OO 

17 The metric at the w — > limit is transformed to ds 2 = —r 2 d,T 2 + r 2 (T 2 dy 2 + dx\) + dr 2 /r 2 by the 
coordinate transformation t = t + -\- \/r. Further boost transformation in the (r, y) directions makes the 
metric the standard Schwarzschild-type pure AdS metric. 
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4.2 Late-time approximation 



We need to justify the geometry (J3TJ) by showing that it is a solution to the Einstein's 
equation ([9]) within an appropriate approximation. We can again make an analogy with 
Ref. [lj to define the approximation we employ. 

We should employ an approximation in which the expansion rate of the fluid is slow 
enough and we can use hydrodynamics. The expansion rate of the Bjorken fluid becomes 
slower and slower along the time evolution; this means that we need to take a large-r 
limit. What the authors of Ref. jl] has found is that we should take the large-r limit 
with v = z/t 1 / 3 kept fixed. One observation is that the naive location of the horizon 
becomes a constant v — Vo in ( |28l) at the leading order if we use (r, t>)-coordinates instead 
of (r, ^-coordinates. The expansion parameter was found to be r~ 2//3 in Ref. [5] by taking 
the viscous effect into account. 

Let us follow the same philosophy to define the late-time approximation on Eddington- 
Finkelstein type coordinates. We introduce a new coordinate variable u which is defined 
by 

u = rrl /S (32) 



so that the naive location of the horizon becomes u =const£j We also define the late-time 
expansion as an expansion with respect to r + with u kept fixed. 

4.3 Summary of our proposal 

We summarize the above discussions, and propose a procedure to construct the dual ge- 
ometry in the late-time regime on the Eddington-Finkelstein type coordinates. 

We propose the following parametrization of the 5d metric: 

ds 2 = -r 2 adrl + 2dr + dr + r 2 T 2 + e 2h - 2c 1 1 + — ^ J dy 2 + r 2 e c dx 2 ± , (33) 

where we have used u defined at Eq. (|32|) to make the order counting transparent. Notice 
that g T+ r + is not parametrized in an exponential form. The parameters a, b, c are expanded 
as follows: 

a(r + , u) = a (u) + ai(u)r^ 2/3 + a 2 (u)r^ 4 ' 3 + a 3 (u)r~ 2 + 0(r+ 8/3 ), (34) 

6(r+,u) = b (u) + hiu)^ 3 + b 2 (u)r- 4/3 + h{u)r~ 2 + 0(r; 8/3 ), (35) 

c(t + ,u) = c (u) + Ci(u)t+ 2/3 + c 2 (u)t^ /3 + c 3 (u)t~ 2 + 0(r+ 8/3 ), (36) 



18 The "naive location of the horizon" means the position where the (t + ,t + ) component of the metric 
across zero. 
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1/3 

where u = rr^ is kept fixed. We solve the 5d Einstein's equation order by order in the 
large- r+ expansion to determine a( n ),6( n ),C( n ). 

The boundary condition we have mentioned around ( TT2|) are equivalent to 

<^|«=oo J-j "|u=oo c| M=00 u. yo I ) 

The stress tensor is identified with ( JT5I) . f[T6l) and (TT7|) by the methods we have presented 
in Section 13.11 and in Section 13.21 

Going through the above procedure, we can explicitly show that (13T1) is a solution to the 
5d Einstein's equation at the leading order of the late-time approximation whose boundary 
condition match the Bjorken flow. We shall show how it works explicitly in Section \5\ and 
in Section [HI Now a few comments are in order: 

• One may be tempted to define g T+T+ in an exponential form like g T+T+ = —r 2 e a . 
However, this is not an appropriate parametrization since g T+T+ cannot be positive 
as far as a is real, despite the fact that g T+T+ must be positive inside the horizon. In 
other words, the late-time expansion of a fails around the horizon. 

• It is quite natural to define g yy and g x±x± by using the exponential forms e 2b ~ 2c and 
e c , because g yy and g x± x± have to be always positive. To see this, suppose that g yy 
reaches zero at a certain value of r in the bulk, for example. Then, the y direction 
shrinks to a point, and different points on the boundary (with the same values of r+, 
x 2 , x 3 but not for y) are mapped to a single point there. The map between the bulk 
and the boundary is ill-defined in this case. The same logic works for g x±x± . 

• We have not yet fixed all the gauge degree of freedom at the metric (l33l . One 
finds that the off-diagonal component 2dr + dr is maintained under the coordinate 
transformation: 

r^r + f(r + ), (38) 

where /(t+) is a function of r+. We shall see explicitly in the next section that the 
un-fixed gauge degree of freedom comes into the solution as an un-fixed integration 
constant. We can use the un-fixed gauge degree of freedom for consistency check; we 
shall find that all the physical quantities are independent of the gauge choice. 

• We can introduce b = 6+log[l + l/('ur 2//3 )] and parametrize g yy = r 2 r + e 2b ~ 2c ; we could 
have started by using b and determine it order by order. However, the advantage 
of our parametrization is that a part of b is already re-summed to all orders in the 
late-time expansions in the form of log[l + l/(ur 2 / 3 )], so that the reduction to exact 
pure AdS is manifest for the empty fluid. (See also the discussion in Section H~T1 ) 
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5 The late-time geometry 

We construct and analyze the dual geometry based on our proposal to the second order of 
the late-time expansion. 

5.1 Zeroth order 

The solution to the Einstein's equation at the zeroth order (leading order) of the late-time 
approximation is given by 

(1 - &/«) 4 - w 4 u~ 4 
a {u) = — 



(1 - W«) 2 ' 
b (u) =31og(l-£ /u), 
c d (w) = 21og(l-£o/u). (39) 

Here, we have already fixed some integration constants so that the geometry matches our 
boundary conditions. Notice that the contribution of l/(w + ) in g yy in ( J33l) has to be 
ignored at this order. £o is an integration constant which cannot be fixed by the boundary 
conditions: £o is a remaining gauge degree of freedom. Indeed, the contribution of £o is 
absorbed by the coordinate transformation 

u^u + ^ + O(r- 2 ^). (40) 

The solution (1391) reproduces the correct boundary metric and the stress tensor. We 
exhibit explicitly the stress tensor of the fluid that is read from the metric at the leading 
order: 

4 i i 4 

T T+T+ =E^, T yy = -Ew'rf , T xx = -E-^, (41) 



where E is defined in ( 1271) . Let us define 

e = Ew\ (42) 

then the stress tensor completely matches (fl~5l) . (TT6l) and (TTTJ) to the leading order. The 
physical meaning of the free parameter w is that it determines the overall factor of the 
energy density. 

As a consistency check, let us compute the Kretschmann scalar. We obtain 

(R„ P x) 2 = 8 (5 + -t^te) + °( r + 2/3 )- ( 43 ) 
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Now we choose the gauge degree of freedom in such a way that the singularity of the 
Kretschmann scalar is located at the origin; we choose £o — 0- Then the final solution at 
the zeroth order is 

ao(u) = 1 — w 4 u~ 4 , 
b (u) = Co(u) = 0, (44) 



which is manifestly regular except at the origin. This agrees with the metric (131]) we have 
anticipated. To make everything consistent, we need an event horizon which covers the 
physical singularity at the origin. We shall discuss this problem in detail in Section [71 

5.2 First order 

The first-order (the sub-leading order) solution is given by 

o M i C "\„,4 i c „,,4 o^, „.„,A 



a\[u) = 


3 




u° 




h(u) = 


6 + 1 

u 






ci(u) = 


- I dx— 

3 Voo x 


x VO i / 1 4 -4\ 

log 1 — W u ) — 

-w 4 2 5V ; 


26 
3u 




1 

3w 


arctar 


... 7T 1 ( U -W\ 






i[U/W) 1 lOg 1 

2 2 \u + w J 






m 

' 2 


log(l 


4 -4\ 26 
-W U )--r-, 

6u 





(45) 

where 6 i s an integration constant which is not fixed by the boundary data. We can 
show that 6 is again a gauge degree of freedom which can be absorbed by the following 
coordinate transformation: 

----J^ + 0(T-n (46) 

Notice that 6 gives the first-order contribution to the transformation in the late-time 
expansion. One useful gauge choice is 6 — — 1- Then a l5 b\ and c\ go to zero at the limit 
of w — > 0, and the geometry is manifestly reduced to pure AdS at w — 00 

Let us check the regularity of the geometry. The Kretschmann scalar to the first order 

is 

(Lf \2 *(k i 9nj8 \ i 96tf7 8 (2^ - 3r/ M) 2/3 nf -4fc (A7 , 



19 However, this does not give any constraint for the gauge. If we choose another gauge, the reduction 
to AdS is still realized order by order in the late-time approximation. 
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and (R^ vp \) 2 is singular only at the origin. 

Note that C\(u) is singular at u = w in general but there is a unique choice 

77o = l/(3w), (48) 

which makes the metric regular except at the origin. Indeed, r/ = l/(3w) is requested 
by the regularity of the geometry in the following way. Let us consider a Riemann tensor 
projected onto an orthonormal basis. One useful component of the projected Riemann 
tensor is R y N^N U where y denotes the rapidity direction and we take the sum only over 
(jl, v. N 11 is a space-like unit vector: 

N" = ~(l,0, OA 1 -^)- (49) 

N^ forms an orthonormal basis together with a time-like vector T M = —4= I — 1, 0, 0, 0, ~ r 2 a+2 j 
on the (r + ,r) plane. We find thato 

W"" = liv^ ("° - WWj) + 0(T ~ n (50) 

and this component is singular at u = w unless rjo = l/(3w). Now, our vectors N^ and T M 
are regular at the vicinity of u = w hence all the components of the projected Riemann 
tensor need to be finite in order to realize a regular geometry. Some readers may wonder 
why the projected Riemann tensor can judge the regularity even though it is not a scalar. 
We provide a detailed explanation in Appendix [Ql 

To conclude, we have shown that, at the first order, the regularity of the dual geometry 
at u = w determines i] to be l/(3w) uniquely. In fact, i] = l/(3iu) corresponds to the 
famous result rj/s = 1/(4:%) [TO] where s is the entropy density. (See Appendix [El) In 
the previous work [10], the condition (148]) was obtained from the condition that (R^upx) 2 
be regular at the second order which is next to ours. (See next subsection.) The reason 
why they have not see the singularity in (R^px) 2 at the first order is due to a non-trivial 
cancellation among the components of the Riemann tensor. 



20 The first-order contribution to Ry, nilJ N^N v is at the order of r°. 
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5.3 Second order 

The second-order solution is given by 

(u 4 - 3w 4 ) H 4 (u 3 - 3w A r] ) a 2 {u 4 + w 4 ) & 



a 2 {u) 



9u 6 9w> 3m 5 

(u 4 - 2w 3 u + w 4 ) (9w 2 r] 2 - 1) 



12u 5 w 



log(u — w) 



(u 4 + 2w 3 u + w 4 ) (9w 2 ti 2 - 1) , . 

(u 4 + w 4 ) (9w 2 ti 2 + 1) /u 

arctan — 



6m 5 w Vu? 

+ -L log (u 2 + w 2 ) 



3 A 



3?7o (3«(12 logw + 5)?7o + 4) w 4 + 4 (3wAw 4 + u 

18u 5 

2u A bu z u A \ u w 

(SwrjQ — 1) (2u — 3w + 3(4w — 3w)wr] ) 
+ ^ l0g(M " W) 

(3wi]o + 1) (— 2u — 3w + 3w(4u + 3w)r] ) 
+ 2W log( " + W) 

+ 1 fiat + J.) i og ( u . + „») + 9ro2 "o- 2 ""° + 1 arct m 

(6(w 4 - 5u 4 ) r] w 4 + 4w 3 (m 4 + w 4 )) fr 2^£ 26 
C2(MJ " 9(^-^4)2 + 9 M 3 + 3« 2 

?7o (12wr/oM 5 — 6wu 4 + n {u 4 — w 4 ) u + 2w 5 ) w 3 

3 (u 5 — uw 4 ) 
Arinu 2 log-u 3r7 -u 3 + w 2 , , , 9 x 7ru 3 — 3u> (4Au> 4 + u 

+ —T-* ^77 T I ^g ( W + IU ) 

q f„,4 „„^ n„,5 n„,„„4 to V / 



(51) 



(52) 



3(w 4 -w 4 ) 9u 5 -9uw 4 v ' 9(u 5 -uw 4 )w 

(3wrj — 1) ((u + w) {u 2 — 2wu + 3w 2 ) — 9(w — w)w {u 2 + w 2 ) i] Q ) 

36u 2 (u — w) {u 2 + w 2 ) w 
(3wi] + 1) ((u - w) (u 2 + 2wu + 3w 2 ) + 9w(u + w) (u 2 + w 2 ) rj ) 

36u 2 (u + w) (u 2 + u> 2 ) w 

u 4 + 3u> 4 — 3w 2 r]o (4uw 2 + 9 (u 4 — w 4 ) r?o) /u\ , . 

arctan I — I , (53) 



log(w — w) 
log(w + w) 



18u 2 (u 4 — w 4 )w \w 

where ^2 is a new integration constant which is a gauge degree of freedom at the second 
order. We can absorb the contribution of £2 by the following coordinate transformation: 

^u-^ + 0(r~ 2 ). (54) 
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A = (A? — Vo T n) IS a combination of the second-order transport coefficients. 02(11) is too 
complicated to present here, and we have exhibited c' 2 (u) instead. (The prime denotes 
u-derivative.) An additional integration constant that appears in 02(11) is fixed by the 
boundary condition C2(u)\ u=00 = 0. 

The second-order contribution to the Kretschmann scalar, R%) T + > can be expanded 
around u = w in the following way: 

^ = l^-^) _ 8^-1) + 

v ' 3[u — w) 1 3{u — w) 

The regularity of the Kretschmann scalar requests 7/0 = l/(3w). Notice that the singulari- 
ties in a 2 and b 2 disappear at r/ = l/(3w). 

Note again that, if we set r] = l/(3w), the coefficient c' 2 is expanded around u = w as 

l + 6^A-log2 

2 18u> 2 (u-w) v ; v ; 

This means that the potential singularity at u = w in c' 2 (hence in C2) disappears if we 
set A = ~ 1 6 l w ° g2 together with r] = l/(3w). Indeed, this value of A is requested by the 
regularity of the geometry as follows. We find that 



R\ yv N»N» 



1 + 6w 2 A - log 2 1 + 6u; 2 A - log 2 
3Qw 2 (u — w) 2 l%w z (u — w) 



+ 0((u-wf) r- 2/6 + 0(r^ /s ) } (57) 

after substituting r] = l/(3w). Therefore, we need 

_ -l + log2 

A = A = -^-o > I 58 ) 

for the regularity of the geometry at u — w. The same condition can be obtained from the 
regularity of the Kretschmann scalar at the third order, which is next to ours [8j [HI [15] . 
We present the details in Appendix [0 



5.4 Summary of the present section 

It is better to summarize what we have found in the present section, before starting more 
general analysis in the next section. We have found the following facts to the second order 
of the late-time expansion: 

• a,b, c and their arbitrary-order derivatives are regular except at the origin if we choose 
the transport coefficientcj appropriately. Although we have not demonstrated the 
regularity of the derivatives, one can explicitly check their regularity as well. 



21 More precisely, the combination of the transport coefficients which appears in the stress tensor. 
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Actually, the foregoing choice of the transport coefficients is a sufficient condition for 
the regularity of the geometry at u 7^ 0. One finds that the inverse metric and their 
arbitrary-order derivatives are also regular at u 7^ if a, b, c and their arbitrary-order 
derivatives are regular. This is due to the nature of the Eddington-Finkelstein type 
metric. If the metric, the inverse metric, and their arbitrary-order derivatives are 
regular, all the curvature invariants are regular. 

One may worry that the metric is divergent at the boundary because of the presence 
of the factor r 2 even though a, b, c are regular there. However, we can explicitly 
show that the expansions of a, b, c around the boundary are 1/u expansions which 
start at the order of 1/u (or higher). Therefore, the geometry at the vicinity of the 
boundary is always AdS. The divergence due to the r 2 factor is just what we have in 
the pure AdS geometry and it is harmless. 

We have found that the above choice of the transport coefficients is also a necessary 
condition to have a regular geometry at u ^ 0. If we take another value of the 
transport coefficient, the projected Riemann tensor becomes singular. 

As a conclusion, the regularity of the dual geometry except at the origin determines 
(the combination of) the transport coefficients uniquely. 



6 Regularity of dual geometry for all orders 

In this section, we generalize the conclusion of the previous section to all orders. We show 
that 

1) We can make the dual geometry regular except at the origin by choosing the stress 
tensor (the combination of the transport coefficients) appropriately, at the arbitrary 
order in the late-time expansion. 

2) The choice of the stress tensor is also a necessary condition for the regularity of the 
geometry at the given order. If we take another value for the stress tensor, the geometry 
has another singularity in addition to that at the origin. 

To show 1) above, it is sufficient to show the regularity of a n , b n , c n and their arbitrary- 
order derivatives for all n, as we have discussed in Section 15.41 If a n , b n , c n and their u- 
derivatives are regular for all n, it is obvious that r + -derivatives of the metric never create 
singularity; what we need to show is the regularity of a n , b n , c n and their w-derivatives 
for all n. We may use "derivative" as the meaning of "^-derivative" below, if it is not 
confusing. For simplicity, we may also use a term "regular/regularity" as the meaning of 
"regular/regularity at «^ 0" in this section. 

We use induction for the proof. The outline is the following. We begin with the 
assumption that a^, bk, Ck and their arbitrary-order derivatives are regular for k < n. We 
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also assume that the expansions of a*,, bk, Ck around the boundary start at the order of 1/u 
or less singular order. Then, the Einstein's equation tells us that b' n and 6" are regular. We 
can generalize the statement to the regularity of b n and its arbitrary-order derivatives by 
integrating or differentiating the equation. We can also prove the regularity of a n and its 
arbitrary-order derivatives in a similar way, by using the Einstein's equation. The proof for 
c n is more complicated since we encounter a potential singularity. However, we find that 
it is always possible to remove the singularity by an appropriate choice of the integration 
constant in a n , which corresponds to the n-th order contribution to the stress tensor. This 
matches our experience; the new transport coefficients 770 and A have been determined by 
requesting the regularity of c\ and C2, respectively. Since we have already shown that our 
starting assumption is valid to the second order, the regularity (under the appropriate 
choice of the transport coefficients) for all order is proved by induction. The proof of the 
statement 2) above shall be given by using the regularity condition for c n . 

In this section, we introduce f = r + and we switch to (f , u) coordinates from the 
{r + ,r) coordinates. Now the late-time expansion is the expansion with respect to f. The 
relationship between the two coordinate systems are summarized in Appendix [Hi We 
define 

a = (l-w 4 u- 4 )+A(f,u), 

b = B{f,u), 

c = C(f,u), (59) 

where A, B, and C contain the all-order contributions starting at the order of f . (Recall 
that 6 = Co = 0.) We write -^-^A as A^ hj \ and -j^A as A' for simplicity. 

6.1 Regularity of b n 

We begin with b n (u). The (f, u) component of the Einstein's equatioio is given by 
( w 2 S /y = -1 (2u{B^f + 3u(C^) 2 - iuC^B^A 

,(0,1)^2 _ 4^(0,1)5(0,1) + 3 (c(o,i))2 + 4 M -i C (o,i) + 2jB (o, 2 A (6Q) 



-M2( B < 



2 

The n-th order contribution at the left-hand side is (u 2 b' n )'f n . One finds that the right- 
hand side at the same order is given by using only bk, Ck with k < n and their derivatives, 
hence regular at u 7^ by assumption. Then we conclude that b n , b' n and b" n are regular 
at m ^ since the integration of the right-hand side over u has no chance to create a 
singularity] 23 ! The regularity at the boundary is confirmed if one counts the power of u 



22 More precisely, the equation coming from the (f, u) component of the Einstein tensor where the 
first component (f) is raised and the second one (u) is lowered. We follow the same notation for other 
components of the Einstein's equation. 

23 Notice that our boundary condition is b n (u)\ u=oc = hence we do not impose a singular boundary 
condition. 
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by taking account of the fact that bk, Ck with k < n are 0(l/u) at the boundary. More 
explicitly, one finds that the right-hand side of (I60j) is 0(l/u 2 ) at the boundary. By 
integrating ( |60l) . we can immediately conclude that the 1/u expansion of b n starts at the 
order of 1/u, and the coefficient of the 1/n-term is an integration constant as we have seen 
in 6 , b\ and 62- We have shown the regularity of b n , b' n and b" n so far. We can iterate 
the above discussion by differentiating (|60|) with respect to u to show the regularity of the 
arbitrary-order derivatives of b n (u). 



6.2 Regularity of a n (u) 

The regularity of a n (u) is shown almost in a parallel way with what we did for b n (u). The 
(f, f) component of the Einstein's equation is given by 

(^)'=T^n, (6i) 

where the explicit representation of F T f is given in Appendix iGl 

The n-th order contribution at the left-hand side is (u 4 a n )'f n . We find that the contri- 
bution of F T f at the same order is given by using only b' n ; ctk, bk, Ck with k < n; and their 
derivatives. Since the regularity of b' n and its derivatives are already shown, the n-th order 
contribution at the right-hand side is regular. Therefore, we can conclude that a n , a' n are 
regular at u 7^ 0. The regularity at the boundary is confirmed in the following way. One 
finds that the right-hand side of (loTj) is 0{u 2 ) at the boundary, just by counting the power 
of u. Then integration of (ISTj) tells us that a n is 0(l/u) at the boundary. We can repeat 
the analysis by differentiating (16TJ) with respect to u to reach the conclusion that a n and 
its arbitrary-order derivatives are regular at u 7^ 0. 



6.3 Regularity of c n (u) 

We obtain the following equation from the (u, f ) component of the Einstein's equation: 

3u 2 {An 4 + u'- w') {2B^ - 3C^) = tf + f 2 f 2 + f 3 f* + f 4 f\ (62) 

where the explicit forms of /1, f'2, fs and f'4 are given in Appendix [G] The n-th order 
contribution to the left-hand side is 

n-l 

3u 2 (« 4 - w 4 ) (2b' n - 3<J f" + 3u e f » Y, o* (2&U - K- k ) ■ (63) 

fc=i 

We can easily see, by counting the number of f derivatives, that the n-th order contribution 
from f 2 , fa, fi contains only au, bk, Ck with k < n, and their derivatives; their contributions 
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are regular except at the origin and the boundary. The contribution from /i to the n-th 
order is 

f x ? = -3nuf n \ - 3u 4 a n - 2u(u* - w 4 )b' n + Aw% n 

+ (regular terms) f", (64) 

where (regular terms) denotes the terms which contain only Ofe, bk, Ck with k < n and their 
derivatives. Combining the above results, we obtain 



2 n 
■~(n-l)b' n + 

3 3n(u 4 — wrj 



c 'n = -o ( n - l ) h 'n + 1T7—A 7T - 3m a n + 4w b n + f TOg (u) , (65) 



where / reg and its arbitrary-order derivatives are regular except at the origin and the 
boundary. We have already shown the regularity of b' n {u). 

The regularity of the second term at the boundary can be explicitly confirmed. We 
should divide /i, fi-, jz-, f± and the second term of (1631) by u 2 (m 4 — w A ) and count the power 
of u\ they are 0(l/u 2 ) at the boundary. Then the right-hand side of (165]) is 0(l/u 2 ) and we 
conclude that c n at the boundary is 0(l/u). However, the second term in (1651) is potentially 
divergent at u — w. From our experience, we expect that we need to choose a new 
integration constant appropriately to make c n regular. The condition for the integration 
constant (that must be related to a combination of the n-th order transport coefficients) 
is given by the regularity of the second term of (|65|) . Namely, the expansion of — 3w 4 a n + 
4w 4 b n + / reg around u = w has to start at the order of (u—w) or higher. From the regularity 
of a n , b n , / reg and their arbitrary-order derivatives, we can write 

a n (u) = C an +0(u-w), 
b n (u) = C bn + 0(u-w), 
/reg(li) = C vcg + 0(u-w), (66) 

where C an , Cb n , C rcg are constants. Then the condition for the regularity is 

- 3w A C an + Aw 4 C bn + C rcg = 0. (67) 

The point is that C an is determined by the boundary metric and the stress tensor. To 
see this, let us go back to (IBTl) and consider the integration constant in a n . The n-th order 
contribution to (16TT) is a linear differential equation of a n and the integration constant 
comes only as a coefficient in the complementary function of the homogeneous equation. 
For our case, the complementary function is 



a (4) 

O) = -=r, (6* 



where a n is the integration constant. One should notice that the integration constant is 
identified with the n-th order coefficient of T TT through 

«i 4) = -4"V, (69) 
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as we have seen in Section [3T2l (Notice that a^ = Y2k a k r + •) 

In the solution to the inhomogeneous equation, we may have other 0(u~ 4 )-contributions 
that originate from the right-hand side of (l6Tf . However, they are independent of the n- 
th order integration constant: they do not carry any information on €q , . Therefore we 
conclude that a n contains e^ in a linear form in the coefficient of the 0(w~ 4 )-term. This 
guarantees that the 0(l)-part of the expansion of a n around u = w carries eg i n such a 
way that 

a n (u) = —Co + constant + 0(u — w). (70) 

Namely, the information of T TT "propagates" from the boundary to u = w in a linear way 
without vanishing. Therefore, C a „ contains a term which is proportional to e , we can 
always adjust e so that the regularity condition fl67j) holds. The linear dependence on e^ 
is explicitly seen, for example, in ( |56j) for c' 2 . Once the regularity of the right-hand side of 
( 165]) is achieved, it is straightforward to show the regularity of c n (u) and its arbitrary-order 
derivatives at u 7^ 0. 

Of course, the regularity condition (!67l) does not depend on the gauge choice; once 
we achieve the regularity at a particular gauge choice, the regularity of any curvature 
invariants does not affected by the coordinate transformation. As a consistency check, we 
can explicitly see the invariance of the regularity condition under the n-th order gauge 
transformation 

Since / reg contains only the lower-order contributions, C reg is invariant under the n-th order 
transformation. The invariance of — 3w 4 C an + 4w A Cb n is shown in the following way. The 
transformation (1TT1) induces 

a n > a n -\- , o n > o n -\- , c n > c n . { 1 z.) 

6u° u 6u z 

Then, — 3u 4 a n + Aw 4 b n is transformed to 

- Zu 4 a n + AwX + ^ ^ n . (73) 

u 

Therefore the regularity condition is invariant under (171]) . 

6.4 Regularity of c' n as necessary condition 

What we have shown so far is the statement 1) we have presented at the beginning of this 
section. The regularity of c' k for k < n so far is a sufficient condition for the regularity of 
the n-the order geometry. Here, we show that the regularity of c' n is indeed a necessary 
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condition for the regularity of the geometry; we show that the n-th order geometry is 
singular if c' n is singular. The quantity we shall examine is the Riemann tensor projected 
on the orthonormal bases: RP tiyv N ,J 'N v . For the metric given in (l59l) . we obtain 

R\,^= 12 J (f + u) t>>^ (74) 

where we have used the same unit vector as (|49l) . hi are functions of u which are explicitly 
given in Appendix [Gj Let us assume for k < n that the regularity of dk, bk, Ck and their 
derivatives have been already achieved by choosing the k-th order transport coefficients 
appropriately. Then we have shown that a n , b n and their derivatives are also regular. Let 
us take only the potentially divergent contribution of the n-th order metric from (TM]) r 4 l 
We find that hi with i > 1 have only regular contributions while ho contains the following 
potentially divergent part: 

d i + \C (75) 

We can check that ( 175]) for n = 1 and n = 2 reproduce the correct divergent pieces of ( !50l) 
and ( 1571) . respectively. 

Therefore, we have shown that the regularity condition ( 1671) is a necessary condition 
for the regularity of R v N^N U at the n-th order; hence (l67j) is a necessary condition for 
the regularity of the n-th order geometry. Combining all the analysis in this section, we 
conclude that the regularity of the geometry at the n-th order uniquely determines e 
which is the n-th order component of T TT . Furthermore, such an appropriate choice of e 
exists for all n. 



7 Apparent horizons 

We have obtained the late-time geometry explicitly up to the third order and we have 
found that the regularity at u ^ is achieved by choosing the correct transport coefficients 
of the fluid. We have also shown that the geometry can be made regular with appropriate 
choice of the transport coefficients to arbitrary higher order. However, this is not enough 
to show that the dual geometry is healthy: we need to show the presence of the event 
horizon which covers the physical singularity at the origin 25 !. 

We show the presence of an apparent horizon instead of that of the event horizon, 
because an examination of the existence of the event horizon in a time- dependent setup is 
not easy. The presence of the apparent horizon is a sufficient condition for the presence of 
the event horizon hence we can prove the absence of a naked singularity [20] . 



24 Recall that the contribution of the n-th order metric to R v fJ/yu N> 1 N L ' is (^(f™^ 1 ). 
25 The physical singularity at the origin could not be seen in the previous works based on the FG 
coordinates, since the coordinates do not cover the region around the origin. See also Section [2j 
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7.1 Definition of apparent horizon 

We define the apparent horizon based on the double-null formalism (36]. (See also Ref. [37] , 
for example.) We foliate the five-dimensional spacetime by null-hypersurfaces E each of 
which is parameterized by a scalar ^ ± , respectively. Let us consider normal 1-forms to 
E which we define n = — d£ . The 1-forms have the null character: g~ x (n ,n ) = 

g^ n ± n ± = 0. 

The normal 1-forms rt^dx^ on our geometry fl33|) on the (r + ,y, X2, x^,r) coordinates 
are given by 

n; = F-(l,0,0), 

n+ = F + (r 2 a,0,-2), (76) 

where the overall normalizations F ± are determined by the integrability conditions g^g^) = 
0. Next, we define null normal vectors (l±) to E^ 1 by l± = e~* g~ 1 {n T ) which are given 
explicitly by 

ft 1 

2F~ 2F~ y ' 

'" K L -(0,0,1), (77) 



2F+ 2F+ 



where we have defined 



e / = -g-\n + ,n-) = -g^n+n' = 2F + F~. (78) 

We can easily check that g(l+, 1-) = —e~f and g(l±, l±) = 0. 
The (null normal) expansions 6± are defined by 

6»± = £±bg/x, (79) 

where C± are the Lie derivatives along l±. Here, /1 is the volume element of the intersec- 
tion!^ of the null hyper surfaces: 

/1 = r^T+e\ (80) 

where b = 6+log ( 1 + — ) . The most important quantity we need to define is G = e^9 + 9_. 
Since F + F~ in e* cancels with (F + F^)^ 1 in 9+9-, is simply given by \9 + 9- where 

0± = 4tlogj/, (81) 



26 The intersection for the present case is the 3d surface where t+ and r are constants, and it is spanned 
by y, x 2 , x 3 . 
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and L± are the Lie derivatives along l±. This means that we do not need to determine F ± 
explicitly to compute 0. 

The trapped region is where > and the un-trapped region is < 0. The apparent 
horizon is the boundary of the two regions; the location of the apparent horizon is given 
by solving 

= 0. (82) 

We do not define the location of the apparent horizon merely by 9 + = or #_ = 0, since 
6± are not invariant under relabellings of the scalars C,± i— > C±(£±) while is invariant [38] . 



7.2 Apparent horizon on the late-time geometry 

Let us compute the location of the apparent horizon (if it exists) on our dual geometry 
order by order in the late-time expansion. for our geometry is expanded with respect to 

-2/3 

r + : 

= O + ir ; 2/3 + 2 r; 4/3 + 0(r- 2 ). (83) 

The location of the apparent horizon (uh) at the leading order is given by solving ©o = 0. 
The location to the first-order is given from 0o + ©iT + = 0. Then it is consistent to 
expand the position of the apparent horizon with respective to r + : 

u H = u + «ir; 2/3 + u 2 r~ 4/3 + 0{t- 2 ). (84) 

We determine uh order by order. 

7.2.1 Zeroth order 

We find 

0o = -~(l-u o -V), (85) 



and we obtain 



uq = w. (86) 



We can easily see the presence of the trapped region since ©o is positive if u$ < w; u$ = w 
is indeed the boundary of the trapped region and the un-trapped region. 

Notice that it is technically important that our metric is regular at u — w. If the 
metric were singular there, the trapped region and the un-trapped region are not described 
by a single coordinate patch. This prevents us from rigorous proof of the presence of the 
apparent horizon unless we find a better coordinate on which we can show that the two 
regions are really smoothly connected. 
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7.2.2 First order and second order 

Let us substitute u$ = w and compute u\. Now = Q\T + + 0(t + ' ) since the zeroth- 
order contribution vanishes by virtue of uq = w. We find 

1 = (6wi + %tw-2^i), (87) 

w 

and we obtain 

6 



^o , 6 


1 


— w H = 


i — 


2 3 


6 



, 3 



where we put the regularity condition ( Hall in the last step. The contribution of £i can be 
absorbed by the coordinate transformation ( 1461) . of course. 

Let us proceed to the second order. We obtain 

12Xw 2 - 72u 2 w + 24&w - 6 log 2 - 3vr + 10 
® 2 = 4^ ' (89) 

after substituting (188|) . Then we find 

12Aw 2 -61og2-37r + 10 6 

u 2 = 1 

72w 3 

= 8-37r-41og2 fr 

72^; + 3 ' l j 

where we put the regularity condition (1381 in the last step. Again, the contribution of 
£2 represents the degree of freedom of the coordinate transformation, and it is absorbed 
by(j5j 



In the above computations, we did not encounter any difficulty like we have pointed out 
in Appendix [A] for FG coordinates; we can compute the location of the apparent horizon 
in a systematic way. The presence of the trapped region is also very clear. The above 
results show that we do have an apparent horizon (hence an event horizon) which covers 
the physical singularity at the origin. We have proved that the singularity at the origin is 
not a naked singularity hence our dual geometry is totally healthy. Furthermore, we have 
shown that the dual geometry is really a dynamical black hole. (The non-staticity of the 
local geometry shall be shown in Section [HI ) 



7.3 Geometrical and hydrodynamical entropy 

Let us compute the volume element (which we denote Ah) of the apparent horizon. We 
obtain 

A H = fi\ u =u H =w 3 - ^r; 2/3 + ^(4 + vr + 12w 2 X + 41og2)r; 4/3 + 0(r- 2 ), (91) 
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where we have already substituted r] = l/(3w), and A is understood to be the physical 
value A . Notice that A H is independent of the gauge choice] 27 ! When we neglect the second 
order term, Eq. (f9T]) implies that the apparent horizon has a smaller area than the r — > oo 
limit value. 

The event horizon coincides with the apparent horizon when the system is independent 
of time. This suggests that (191]) has to agree with the volume element of the event horizon 
at the infinitely far future. Then the leading order contribution in (I9TT) divided by 47rG 5 , 
where G5 is the 5d Newton's constant, must be the entropy density at the infinitely far 
future. In our convention, (A^G^)" 1 is iV 2 /(27r 2 ) and the entropy density per unit rapidity 
at the late-time limit evaluated from fl9Tl) is 



2„,,3 



T4.S 



~2^ 



(92) 



where s is the entropy density per unit physical volume. 



Let us compare the above results with what we obtain from the fluid dynamics. The 
entropy density (per unit rapidity) we get from the hydrodynamics is 



TS 



iV> 3 
2vr 2 



1 _ ^ r -2/s , 3JJJ-2A 4/3 ' 
2 4 



+ 0(t- 



(93) 



The leading-order contribution of (I93p completely agrees with what we expect from the 
volume element of the apparent horizon at the late-time limit (1921) . Indeed, the first-order 
contributions also agree between (J9"T|) and (193]) if we regard r + = r, while it is not the case 
for the second-order contributions. The second-order contribution to the entropy density 
obtained from f]9T]) is larger than that in (l93ll 28 l 

However, the disagreement at the second-order does not necessarily mean a physical 
inconsistency. First, an important fact is that Ah is evaluated at the position of the 
apparent horizon and r + at the horizon can be different from the proper-time at the 
boundary. We have an ambiguity how to map the proper-time at the boundary to the 
horizon. Second, there is still room for discussion whether we can employ the volume 
of the apparent horizon (with an appropriate normalization) as an entropy of the time- 
dependent system or noto Therefore, what we can compare concretely at this stage is 
only the time-independent piece of the entropy density on which we have the complete 
agreement. 



27 If we substitute a wrong value to 770 formally, Ah has a £1 dependence. 

28 Furthermore, the entropy density computed from the area of the event horizon is equal to or more 
larger than that obtained from the apparent horizon, since the event horizon is not located inside the 
apparent horizon. 

29 A formulation of the first law of thermodynamics by using the dynamical apparent horizon for 4d 
geometries is proposed in Ref. [37]. However, its generalization to 5d geometry is not straightforward. 



30 



8 Non-staticity of the local geometry 

In this section, we briefly comment on the non-staticity of our local geometry. Because 
of the time-dependent boundary condition, the geometry we have considered is not glob- 
ally static. However, there are some examples in which a globally non-static geometry is 
actually locally static. One example is the Randall-Sundrum braneworld [39] in a cosmo- 
logical setup where the Friedman-Robertson- Walker universe is realized on a brane. The 
5-dimensional geometry is globally dynamical but locally static in this example [40j. The 
reason is that the geometry has the symmetry of a 3-dimensional constant-curvature space, 
which enables us to apply a generalization of the Birkhoff theorem to ensure the staticity of 
the local geometry. The local staticity is seen in the following way. In the Schwarzschild- 
like coordinate system, the bulk geometry is written in a static form but the brane is 
moving in the radial direction [41, 42J. Actually, the geometry is locally Schwarzschild- 
AdS spacetime [43]. On the other hand, in the Gaussian normal coordinate system, the 
brane position is fixed relative to the coordinate system but the bulk geometry is written 
in a time-dependent form [44, 45j . 

Therefore, it is not a priori clear whether the bulk geometry considered in this paper 
is locally non-static. In the following we shall show that it is indeed locally non-static. 
For this purpose we shall pay attention to evolution of the anisotropy between x and y 
directions. If the symmetry between x and y directions is broken then the generalized 
Birkhoff theorem mentioned above does not apply: we expect the anisotropy to evolve 
unless the boundary condition is very special. (See, for example, Ref. [46J.) 

We obtain the following expansion of components of the Weyl tensor for our geometry: 

°^-- = ^-S^ 2/3+o(t ~ 4/S) - 

These components do not show anisotropy between x and y directions at the limit of 
r —>■ oo. However, for a large but finite r there remains anisotropy if r] ^ 0. This means 
that the anisotropy evolves in time and that the geometry is not locally static under the 
presence of dissipation. (Recall that the dissipation also makes the volume element of the 
apparent horizon to be time dependent.) 



9 Conclusion and discussion 

We have studied a gravity dual of Bjorken fluid at the late-time regime. We point out 
the problems of the Feffer man- Graham coordinates and we propose a recipe to construct a 
dual geometry on Eddington-Finkelstein type coordinates. We have constructed the dual 
geometry explicitly to the second order of the late-time expansion. We have found that 
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the regularity condition uniquely determines the transport coefficients: the shear viscosity 
(148]) and the combination of the second-order transport coefficients (1581) . They agree with 
the results obtained by other methods. 

We have also shown that the regularity of the dual geometry is realized at all orders 
by choosing the transport coefficients appropriately. This means that the logarithmic 
singularity discussed in Ref. [9] is absent from our dual geometry. Our interpretation is 
that the large-r expansion is ill-defined at the vicinity of the singularity (or at the vicinity 
of the "would-be horizon") on the Fefferman-Graham coordinates. We have also proved 
the presence of the apparent horizon (hence the event horizon) on the dual geometry and it 
is shown that the geometry is really a dynamical black hole. The singularity at the origin 
is not a naked singularity and the dual geometry is totally healthy. The metric (with our 
choice of the coordinates) is found to be regular if we use the proper transport coefficients. 
The regularity of the metric was also technically necessary to carry out the analysis of the 
apparent horizon since we need to use a coordinate system which covers both the trapped 
and the un-trapped regions smoothly. 

We can summarize how the hydrodynamics of the 4d YM theory is encoded in the grav- 
ity dual as follows. The hydrodynamics is an effective theory in which the transport coeffi- 
cients are free (un-determined) parameters. We have the hydrodynamic equation, however 
we need the equation of state to solve the hydrodynamic equation, and the equation of 
state is given by the microscopic theory. As is pointed out in Ref. [15], the hydrodynamic 
equation is obtained by the Einstein's equation (around the boundary). The equation of 
state (the traceless condition) is also a consequence of the asymptotically AdS spacetime 
which is ensured by the Einstein's equation and the boundary condition. An interesting 
fact is, on the other hand, that the transport coefficients are determined by the regularity 
around the (apparent) horizon which is deep inside the bulk. 

It is interesting to consider what classifies the properties determined around the bound- 
ary and those determined around the horizon. The traceless property and the conservation 
of the stress tensor hold whether or not the (local) thermal equilibrium is achieved in 
the YM-theory side. However, the concept of the transport coefficients makes sense only 
when the notion of fluid is valid. This tempts us to relate the notion of (local) thermal 
equilibrium with the regularity (or the presence) of the horizon. It is also interesting to 
see how the method to determine the transport coefficients from the regularity is related 
to the Kubo's linear response theory (Kubo formula) and other holographic computations 
(see for example, reviews [II]). We hope that these points will be clarified in the future. 

We have also discussed the proper-time dependence of the entropy density from the 
viewpoint of the dual geometry. It is also interesting to pursue this direction further Jfj For 
example, it is interesting to study a thermodynamic formulation of dynamical black holes 
in asymptotically AdSs geometries by generalizing the work of Ref. [37]. Identification of 
the times at the boundary and the horizon also calls for further consideration. Our model 
provides a consistent setup for the holographic dual of Bjorken flow of M = 4 SYM plasma. 



30 Related works can be found at Refs. [H [H gH H| • 
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The model serves as a concrete well-defined example of time-dependent AdS/CFT, too. We 
hope that the present work sheds some light on the dynamical nature of the time- dependent 
plasma. 

Note added: 

When the present work was at the final stage, we have received a paper [50] which overlaps 
with our results. The first-order solution presented in Ref. [50] corresponds to the gauge 
choice of £1 = in our first-order solution^ Our proposals in the present paper have been 
invented independently. However, we were motivated by Ref. [50] to examine the gauge 
degree of freedom and the regularity of the higher-order geometry. 
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A Failure of late-time approximation on 
Feffer man- Graham coordinates 

We demonstrate how the computation of the location of the apparent horizon fails on the 
FG coordinates. Let us attempt to carry out the computation based on the late-time 
expansion. We assume that the position of the apparent horizon can be expanded in the 



31 Recall that with this gauge choice, it is non-trivial to show the smoothness of the empty limit to the 
AdS space, see discussions in Sections [4.11 and l4~3l 
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following way: 

vh = v + v x t~ 2Iz + v 2 t~ 4/3 + ■■■ . (95) 

Then based on the dual metric obtained in Refs. [H [5j [7J [8] is given by 

= O + lT - 2 / 3 + 2 r" 4 / 3 + • ■ ■ , (96) 



where 



©o = - 



9 (3 ~ vtf 
2(3 + t$" 



_ 108v 3 (3~v^(2v 1 -v fi ) 
(3 + ^ 4 ) 3 

02 = 2(3 + V ( (3 " V ° > (3 + ^ (3 + Cvl ~ V ° + 1UVoV2) 

+ 72(3 - vj})(3 + v 2 )(3 - 5vS)v? - 288v (9 - 12^ + v^ Vl fj 

- 24w 2 (-45 + 45^ + 2^)?7 2 ) . (97) 



Here, C is a constant related to the transport coefficients: 



C = 36(4 2) -H^), (9f 



where r/o is a parameter which characterizes the shear viscosity! 32 ! 

From the expression of 0o, one may conclude that vq = 3 1//4 . However, its justification 
is not clear. The reason is that the approximation in (1971) is not valid around vq = 3 1 / 4 . 
This is due to the following fact: 

• ©o is proportional to (3 — t> 4 ) 2 , 0i is proportional to (3 — t> 4 ) and ©2 has a term 
which does not contain (3 — t> 4 ). This means that an effective expansion parameter 
is (3 — Vo) _1 t~ 2 ' 3 rather than r~ 2 ' 3 (with an appropriate dimensionful coefficient). 

Namely, the effective expansion parameter becomes infinitely large at Vq = 3 1//4 as far as r 
is finite. 

The failure of the approximation is also seen in the following way. If we attempt to 
compute without expanding vu with respect to r~ 2 / 3 , the effective expansion parameter 
of becomes (3 — f 4 I )~ 1 r -2 / 3 . Next, we attempt to substitute (|95|) into the expression of 
0. Then we find that v = 3 1//4 , and 3 — v# is at the order of r~ 2//3 . As a result, the effective 
expansion parameter becomes 0(1) and we are not employing the large-r approximation 
anymore; we need all-order resummation to get a sensible result. 

Therefore, even if we find a region where > by truncating the late-time expansion of 
at some order, we cannot conclude the presence of the trapped region. We can explicitly 
see that the result strongly depends on how we truncate the expansion. 



32 fjo is denoted as 770 in the Refs. 0EHS], and the value 770 = 2 1 / 2 3 3 / 4 corresponds to r\j s — 1/(4?t). 
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B Induced metric and extrinsic curvature 

We define the induced metric 7^ on the regularized boundary (which is an r =const. 
surface) by 

T/w = 9tu> ~ n^nv, (99) 

where g^ is the bulk metric and n M is the outward-pointing unit normal vector to the 
regularized boundary. Notice that f y iai n v = and 7^ is a 5 x 5 matrix. We define 
l^ v — 9 fia 9 u ^'lai3- For our geometry defined in ffTTj) . or in (1331) . the normal vector is 
explicitly given by 



rv= 0,0,0,0,^= , (100) 

where g 00 = r 2 g T+T+ and —r 2 a for (1TTT) and (I3"3"j) . respectively. 
The boundary extrinsic curvature is given by 

K^ = -\{V»n v + V v n») = -^V^V^ + ^Vph a ). (101) 

Here V M is the covariant derivative with respect to 7^, and ^V M is the 5d covariant 
derivative with respect to g^. Notice that K IMU h l/ = 0. We also define K = K^g^ u = 

The boundary Einstein tensor is defined as 

G^ = WRp, - ^ {4) R, (102) 

where the curvature tensors are defined by using 7^ and 7^. They are related to the 5d 
curvature tensors (defined with respect to g^J) through the Gauss equations: 

(4) R«^v = ^RnXpal^l^l^V + K^K^-K^K^ (103) 

MRp, = ^R a ^ ul ^ = l KX l, p lu a{5) R^x, + KK„-K, a Kl, (104) 

(4) i? = (4) ^7^ = {5) R - 2h a n ®R ap + K 2 - K a/3 K a ^ } (105) 

where ( 5 ) represents the 5d quantities and we have put ^ to the curvature tensors associated 
with the induced metric. All (199]) . (110 II) and (11021) are 5x5 matrices. The 4d stress tensor 
is given by using 7^, K^ and G^ in ( ITBl . The stress tensor is defined as a 5 x 5 matrix 
there, but only the 4x4 part has to be taken when we read the 4d stress tensor, of course. 

Another definition of induced metric which may be more familiar to particle theorists 

is 

% = diX^X^, (106) 
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where X^ denote the spacetime coordinates on the regularized boundary and 7^ is defined 
as a 4 x 4 matrix. (Let i, j run from to 3 in this section.) We may write the extrinsic 
curvature as 

K VJ = —diX^X^V^ + (5) V,n M ), (107) 

and K = Kifi l \ where 7 1 - 7 is the inverse of 7^. The stress tensor is also given by using 7^, 
Kij, the boundary Einstein tensor defined with respect to 7^ and the 4x4 counter part 
of (fT3"|) so that everything is written by 4 x 4 matrices. 

Both the above two methods yield the correct 4d stress tensor. In the present paper, 
however, we have used the definition based on the 5x5 matrices rather than the 4x4 
matrices, since it is more convenient for the discussions at Section I3.ll 



C Second-order hydrodynamics of conformal fluid un- 
der Bjorken expansion 



A new second-order hydrodynamics has been proposed by Refs. [HI [15] recently. We review 
it along Ref . [H] . The stress tensor of the fluid can be decomposed aa 33 l 

T^ = eu^u u + PA^ U + IP", (108) 

where u^, e and P are the 4- velocity, the energy density and the pressure of the fluid, 
respectively. The spatial projection A M " is given by A M " = g^ A , + u^u u and IP" is the 
dissipative part. 

The dissipative part in the second-order hydrodynamics proposed in Refs. [T4"l [15] has 
the following expression: 



(D^ + VfV-w) 



IP" = -7?<P" + r]T U 

+ Aicr^ x o v)x + A 2 <x <M A ft" )A + X 3 Q {fM A ft" )A , (109) 



where D = iPV^, 77 is the shear viscosity, Tn is the relaxation time, and Ai, A2, A3 are 
the new second-order transport coefficients introduced in Refs. [TJJ [15] . <t m " is defined as 
<P" = 2^V m m"^ and £P" is the vorticity [H]. The bracket in the indices means 

{ A^ = l - A^ A^(A a/3 + A Pa ) - 1 A^ A a ?A a0 , (110) 

which is the traceless transverse part of the second-rank tensor, projected onto the spatial 
part by A. In Eq. (1109p . we have assumed that the fluid is on a flat spacetime and have 
omitted the curvature dependent part given in Ref. [T4"] . 



33 In this section, /i, v and other indices run from to 3. 
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Let us consider a fluid which undergoes the Bjorken expansion [3j. We set W = 
because it is absent from the Bjorken flow. We choose our coordinate to be the local rest 
frame (r, y,x 2 ,x 3 ) where v? = (1,0,0,0). r and y are the proper-time and the rapidity 
of the fluid respectively, and x 2 , x 3 are the transverse directions to the expansion. The 
dissipative part in this setup is explicitly given by 

/0 \ /0 \ 



W u 



4 r -3 
3' 



V 



■2-r- 1 

3 



1 ) 



+ (v T n - ^i, 



-8 T -4 

9' 



V 



9' 



4 r -2 
9' 



/ 



111) 



and the non-dissipative part is given by diag(e, P/t 2 , P, P). The stress tensor of the con- 
formal fluid is traceless and the equation of state is e = 3P. 

Let us solve the hydrodynamic equation V^T^ = 0, which is explicitly written as 

S(Ai - r]T U ) 



d T 



e + P 4?7 

~~^ + 3^ 



9r 3 



(112) 



From the conformal invariance of the fluid, the proper-time dependence of the transport 
coefficients are given by using that of the energy density: 



V = eo?7o 



i) 



3/4 



Tu 



Tu — 



-1/4 



Ai 



^ U 



1/2 



;ii3) 



where e , r/ , r^ and A° are constants. By using the above equations together with the 
equation of state, the solution to the hydrodynamic equation in the late-time regime is 
obtained to be |14| 



e r 



fo 



where 



-4/3 _ 2r]oT -2 + e (2) r -8/3 + . 



M _ 9t7 2 + 4(A° - Vo r^ 
z ° ~ 6 



;il4) 



(115) 



The higher-order terms denoted by dots are ignored in our approximation. Notice that 
the power —4/3 of r at the leading order is obtained from the first term in (11121) and the 
equation of state (e = 3P). The non-zero components of the energy-momentum tensor in 
the late-time regime are then given by 



T TT /eo 
T yy /e 



-4/3 



2r ]0 r- 2 + e ( i ) r-^ + .. 



■ XJ_X 



./to 



T 
1 

3 

L--V3 - M 2 V-8/3 



.2/3 



-2 Vo +-c, r 
1 



(2)^-2/3 



;ii6) 
;ii7) 

'118) 



3 3 u ' 

where x± denotes x 2 and x 3 . Notice that the indices in the stress tensor have been lowered 
in the foregoing expression. 
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D Riemann tensor on orthonormal basis and singu- 
larity 

Riemann tensors are not coordinate invariant quantities. However, we can conclude that 
a geometry is singular if a component of the Riemann tensor projected onto a (regular) 
orthonormal basis is singular. The projected Riemann tensor is 

where e^ is a vielbein and a, b, c, d range over the five- dimensional (local) Minkowski 
coordinates (which is spanned by the orthonormal basis) o The Riemann tensor is the 
unique lank four covariant tensor out of the metric. Its vielbein component is the resultant 
of the projection onto the local Minkowski spacetime; we still have an ambiguity due to 
the degree of freedom of boost and rotational transformation on the Minkowski spacetime. 
However, the remaining ambiguity does not affect our conclusion in the following sense. 
Suppose that one finds some component of R a bcd is singular at some point in the bulk. The 
singular nature does not change under the remaining boost and rotational transformations 
unless we consider an infinite boost. Therefore, if we find a singular projected Riemann 
tensor^jthe geometry is singular and some (non-trivial) curvature invariant must be singular 
there] 



One useful component for us to see a potential singularity is R y IMyv N^N u as we have 
discussed in Section 15.21 Since our metric is diagonal in the y direction, the vielbein and 
its dual basis cancel each other when they act the upper and the lower indices labeled 
with y. We cannot use the same trick for the components which contains r + and r, since 
they are not diagonal. Therefore, we have considered the projection explicitly by using the 
inner products with vector iV M . Since N^ and T M in Section [531 are regular at u = w, our 
local orthonormal basis is well-defined there. We could have used another basis to reach 
the same conclusion so long as the basis is connected to ours by a finite boost. 



E The ratio of viscosity to entropy-density 

We assume that the static result for the relationship between the energy density and the 
temperature [19] holds in the late-time regime: 

e = -n 2 N^T\ (120) 



34 We have used a, b, c and d for the Minkowski coordinates. 

35 However, we cannot conclude in a opposite way; the geometry may be singular even though all the 
projected components of the Riemann tensor are regular. The reason is that we are considering only to 
the second derivatives of the metric within this discussion, and we cannot judge the curvature invariants 
which contains higher-order derivatives. 
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From the conformality, the entropy density is proportional to N 2 T 3 . The precise coefficient 
is determined by using the first law of thermodynamics dF = —sdT and F = e — Ts, where 
F is the Helmholtz free energy per unit volume. We obtain 

a = yiV c 2 T 3 . (121) 

In our notation, 77 is given by f II 1 3 j) with the obtained value of e at (|42j) . Combining these 
equations, we obtain 

V - = ^-3 Vo w, (122) 

s An 

which yields l/(4*7r) by substituting our result from the requirement of the regularity r] = 
1/(3™). 



F Kretschmann scalar at the third order 

The third-order solutions are too complicated to present here, and we discuss without 
presenting the explicit solutions. For this purpose, we utilize the Einstein equation on the 
(f , u) coordinates. 

We begin with analysis of b' 3 (u) and 0,3(11) around u = w. We find that the analytic 
solution b 3 (u) can be expanded around « = iyin the following way: 

b' 3 (u) = 4 0) + $\u -w) + bf(u - wf + 0((u - wf) 

+ (A - A )log(w - w) fe 0) + B { z ] (u -w) + Bf\u - wf + 0((u - w) 2 )] , (123) 

where 

and 63 s and -B3 s are nonzero constants. Next, let us consider the (r, f) component of 
the Einstein equation. We obtain 

8u 3 o 3 (u) + 7u 4 o' 3 (u) + u b ol(u) + 2(u 4 + w 4 %(u) = f f , (125) 

where 

2^(A-A„) 
3(w — w) 

By combining ( 1123ft . ( 1 125ft and ( 11261) . we conclude that a 3 ', a' 3 and a 3 are less singular than 

(u — w)~ 2 at u = w. 
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Let us show that A = Ao is necessary to achieve the regularity of the Kretschmann 
scalar by using the above observation. We find that the third-order contribution to the 
Kretschmann scalar, -R 2 3 nT+ 2 , is given by 



16w(A-A ) | 16(A-A ) | 



R ®^' 3(„ ,rf 



80o 3 (m) + 40m- 



u — w 
24m; 4 



w 



,, x , 4(u 4 -3w 4 ) „ 

a sW + ~2 fl s( M ) 



ir 






U' 



[u — w 



ir 



-b'iiu). 



(127) 



The full expression is too much complicated and we have expanded the first line around 
u = w. Since 03,03,03,63 and 63 are less singular than (u — w)~ 2 , we have no way to remove 
the singularity at the order of (u — w)~ 2 unless we set A = A . 



G Supplement for the all-order analysis 



F T f in (161]) is given by 



-3m 2 



+f 



A(Au 5 + u 5 -w 4 u)(B^) 2 

+2jB (o,D (A(°- V + IOAm 4 + 10m 4 - 2 (Au 4 + m 4 - m> 4 ) C^u - 6m; 4 ) 

+u (Au A + m 4 - w 4 ) (3{C^) 2 + AB^) ' 

6^(0,1) ^(o,i) + 2 )u 5 - VIA (uB^ (uB<M + 3) + 3) u 4 

-3 (3Au 4 + 5m 4 - 3m; 4 ) {C^fu 2 - 4 (3Au 4 + 5m 4 - 3m; 4 ) B^u 2 
+AB^ (3 (m; 4 - 6m 4 ) + (3mm; 4 - 5m 5 ) B^) 



+a{u (3Au 4 + 5m 4 - 3w 4 ) B™ - 3 (Au 4 + m 4 - w 4 ) }c^ 



1 1 



-12m; 4 



-2fV 



+8fV 



Au(B^) 2 + 3m(C(°' 1 )) 2 + 4C(°' 1 > + AuB^ - 12B™ 
+S (o,i) (q _ 4u(C^ + B™)) - AuB {1 ^ 
(uB^ + 2)B^+uB^ 



'1281 
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The right-hand side of (1621) consists of 

fl = -l2w 4 -3A^(uB^ + l)u 5 

-u \uBW f - 3A^u 4 + 2u 3 + Q (Au 4 + u 4 - w 4 ) (C™ - 5 (1 ' 0) ) 

+3J - 3^'V - 2 (Au 4 + u 4 - w 4 ) B^u + 4wW> 

+C^ [A (2uB™ - 3uC^ -l)u 4 - 3m 4 ) 

+C^ ((u -w){u + w) (u 2 + w 2 ) (2B^ - 3C^) u + w 4 ) }] , (129) 

f 2 = -3A^B^u 5 + 6A^u 4 - QAB^u 4 - AB^u 4 + GAC^u 4 - 6C^u 4 
-C^ ((3Au 4 + 5m 4 - 3w 4 ) (2B™ - 3C^) - 2u 3 ) u 



.£(0,1) (_ 3A (1,<V _ 2u 3 _ 2 ( 3Au l + 5n 4 _ 3w 4) ^(1,0) _ ^(1,0) j j 

-2 (3Am 4 + 5m 4 - 3w 4 ) B^u - 6m? 4 5 (1 ' 0) - 6m? 4 C (1 



u 



o) 



h 



h 



-2u 



2u{B^f + (5 - 2u (B^ - C^ + 2C^)) B™ 

+u (C^ {2B^ - 3C^ + 3C^) - 2B^ + 2B^) 
2u 3 (2(B^) 2 - 4C(i.o) S (i,o) + 3(c(i-o))2 + 33(2,0)) 



(130) 

(131) 
(132) 



The n-th order contributions of F T ~ , f 2 , fa and / 4 contain only b n ; a^, 6fe, C& with k < n; 
and their derivatives. This can be easily seen by counting the number of f derivatives 
and by taking account of the fact that A, B, C are O(f). The power of 1/m around the 
boundary is also readable by taking account that A, B, C are 0(l/u) in the large-u region. 

The right-hand side of (I74p is given by using the following functions: 



he 



h 



-36m 4 



-36m 



tj (£(0,1))" + ( 2 _ 2mC(°< 1 )) B^ - 2C^ 



(133) 
A (o,D ^£(o,i) _ uC ,(o,i) + i) M s + A ^ u ( 5 (o,i)) 2 + ( 4 _ 2uC^) B^ 

-AC^ + u {{C^) 2 + £(°- 2 > - C< ' 2 >) ) + 2}n 4 

+ { - 4C<°'V + (m 4 + m - m; 4 ) (fi^) 2 + (m 4 + m - m; 4 ) (C^) 2 
+2B<°» 1 > (2n 3 - (m 4 + u - w 4 ) C^) + (u 4 + u- w 4 ) (B^ - C^) }m 2 

+2(m 4 + m; 4 )1, (134) 
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24m 4 + (3w 8 - 6w(w 3 + 2)w 4 + m 5 (3w 3 + 20) 
+3Au\Au 4 + 2m 4 + 4m - 2w 4 )\ (B^ 1] ) 2 u 
+ (3w 8 - 6w(w 3 + 2)w 4 + m 5 (3m 3 + 20) + 3Au 4 (Au 4 + 2m 4 + 4m - 2m> 4 )) 

-2 (3w 8 - 6m(m 3 - 2)w 4 + m 5 (3m 3 + 32) + 3Au A (Au A + 2m 4 + 4m - 2m> 4 ) 
+6u e A^)c^ 



+2B^ ( - 3A 2 (uC^ - 1)m 8 + 6A^u 6 + (3m 3 + 32)m 5 
-QA(u A + 2u- w 4 ){uC^ 1] - 1)m 4 - Q(u 3 - 2)w 4 u 
- (3w 8 - 6m(m 3 + 2)m> 4 + u 5 {3u 3 + 20)) C (0>1) « + 3m> 8 



(135) 



3A 2 (C^)\ 8 + 10A (C^ 1 )) V + 7 (C^ 1 )) V + 3A 2 5(°> 2 )m 8 + lOAB^u 8 



+ 7jB (0,2) M 8 _ 3^2^(0,2)^8 _ 1MC7 (0,2) M 8 _ 7^(0,2)^8 _ ^(0,1)^7 



I 8 

. 2 



lGAC^V 



-lGC^V + 8 (C^) V + 8B^u 5 - 8C^u 5 + WC^B^u 5 

_ 16C (0,1) C (1,0) W 5 _ 16S (1,D U 5 + 16C (1,1) W 5 _ q Aw 4 ( C (0,1)) 2 W 4 

-10m; 4 (C^ 1 *) V - 8C^u 4 - 6Aw 4 B^u 4 - IOmAB^m 4 + 6Aw 4 C^u 4 
+10w 4 C^u 4 - 16B (1 ^u 4 + lGC^V + WC<°>V - 16w 4 u 2 
+ (7m 8 + 8m 5 - 10w 4 u 4 + A (3Au 4 + 10u 4 - 6w 4 ) u 4 + 3w 8 ) {B^f 

+3w 8 (C^f + Sm; 8 ^ ' 2 ) - 3w 8 C^ 

-2B^ (3A 2 C^u 8 + 2A ((5m 4 - 3w 4 ) C™ - 4m 3 ) m 4 

-4 (2m 4 + 2 (C^ - B^) u 2 + u- 2w 4 ) u 3 

+ (7m 8 + 8m 5 - 10m; 4 m 4 + 3m; 8 ) C^A] , (136) 
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hi. = 



h- 



-Au' 



3AB^u 5 + AB^u 5 - 3AC^u 5 - AC^u 5 + 3A^B^ 



u' 



(137) 



_ 3A (o,i) c .(i,o) M 5 _ dAB (i,D u 5 _ 6B (i,D u 6 _ 3A^u 4 - 12B^u 2 
-3w 4 B^u + 3w 4 C^u 

+C (0,1) (3 ( A d,0) u 4 + 2 (^4 + ^4 + 2m _ ^ ^(1,0) _ ^(1,0))) _ 4m 3) ^ 
+jB (0,l)| _ 3,4(1-0) M 4 + ^3 + (_ Mm 4 _ gu 4 + g^4) £,(0,1) 

+6 (Au 4 + u 4 + 2u- w 4 ) (C^ - B™) }u 

+Qw 4 B {1 ' 1) u + 6 {Au 4 + u 4 + 2u- w 4 ) C {1 ' l) u 

+ (3Au 5 + Au 5 - 3w 4 u) (B^) 2 + (3Au 5 + 4u 5 - 3w 4 u) (C^) 2 

+12w 4 B^-12w 4 C^], 

h 5 = ^[(BMj^^j^Yii^B^^'-^V + SiMBWti 5 
_ 3A (o,i) c m u s _ 6AB* 1 ' V - lOB^'V + QAB^u 4 

+4jB (l-0) M 4 _ 6AC d,0) u 4 _ 4C -(1,0) U 4 

+C r(o,i) ( 3i4 (i,o) w 4 + 2 (u 3 + (3Am 4 + 5u 4 - 3w 4 ) (B™ - C^))) « 
_ B (o,i) ( 2 C(°>V + SA^V + 2 (m 3 + (3Am 4 + 5m 4 - 3w 4 ) (B™ - C^))) u 
+6w 4 B^u + 2 (3Au 4 + 5u 4 - 3w 4 ) C (1 ^u + 6w 4 B^ - 6w 4 C^°A , (138) 

- 2u (bW>) 2 + (2m (B^ - C^ + 2CW) - 5) bW> + BC™ 
2u (C™ (B^ - C^ + C™) - £ (1,1) + C^ + B< 2 '°) - C< 2 '°>) 1 , (139) 
16m 4 ([B™ - C^f + B^ - C^) . (140) 

One finds that the n-th order contributions from hi with i > 1 are regular by counting the 
number of f derivatives and the power of f in (1741) . 

H Equations in (f , it) coordinates 



->» 



Practically, it is more convenient to work in the following new coordinates to perform the 
f = t, expansion with fixed u = rrl . We list the one forms 
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The Eddington-Finkelstein type metric and its inverse in our ansatz read 

g = - f^a + ^] di 2 + e 2b - 2c (- + lYdy 2 + e c fu 2 dx\ - -dfdu, (142) 

\4r 4 2t 6 J \t / r 

9- 1 = ^ wc j^f d v + S c - x 9A + (r 2 + t9 *■ (l43) 

where da^dx" = ^ M ^^ JV +^ JV ®^ M an d Q M d N = d M®d N +d N ®d M ag ugual 
We take the normal 1-forms n^ 1 as 

n\ l dx M = -F + ((3au 2 + 2«f) df + 4f 2 dw) , 

n M dx M = -F~df. (144) 

We can check that 

g-\n±,n ± )=g NM n%n± d = (145) 

and 

e f = -g-\n + , n-) = -g NM n + N n- M = ^F + F~f A . (146) 

The intrinsic metric h = g + e~f {n + <S>n~ +n~ <S> n + ), or in components Hmn = gMN + 
e ^ { n M n N + n M n ~N)i becomes hfM = h U M = and hij = gij, where i,j, . . . run for 1, 2 
and 3. More explicitly, 

h u = e 2b ~ 2c ( - + 1J , h 22 = h 33 = e c fu 2 , others = 0. (147) 

The perpendicular volume element is 

Vdet {3) h = e b u 2 (u + f). (148) 

The null normal vectors are given as l± = e~J g~ l {n T ), that is, l± = e~$ g MN mJfj 

1 +°m = 7j^A> l - d M = jz{ d f- W + 2f i^ d u) ■ (149) 

The expansion is given as the derivatives l± of the perpendicular volume 

0± * 7S^ (< ± ^5%) = -^^ im + JJfc) «V<« + r)] . (150) 
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I Lie derivatives 

For an unfamiliar reader, we list the basic formulae for the Lie derivatives that is employed 
in Refs. [361 137] in the definition of the expansions. For a 0-form / and for basis of tangent 
and cotangent spaces <9 M and dx 1 *, the corresponding Lie derivatives along a direction 
X = X^d^ are given by, respectively, 

Cxf = Xf = X» (dj) , Cxd^ = - [d^X v ) dy, Cxdx* = {d v X») dx v . (151) 

The Lie derivative for general expressions can be obtained from 

C x (ti ® t 2 ) = (C x ti) ® h + ti ® (£^ 2 ) , (152) 

where ti and £2 are tensor fields of arbitrary types. For a vector F = Y u d u , 

£xY = {X(Y»))d u + Y»(C x dt t ) 

= {X"d li Y v - y<Vn d„ = [X, Y], (153) 

for a 1-form u = u^dx 11 , 

Lxu = (X(u v )) dx v + Up (Cxdx 11 ) 

= {X^u v + u^X") dx\ (154) 

and for a mixed tensor, say, t = t p dx^ ® ct, 

£ x i = (X(y)) do" (8) <9„ + £/ (£*<&') (8) d„ + £/ efcr" (8) (£*#„) 

= (X p d p t^ + t/d^X? - t/d p X u ) dx» ® <9„. (155) 

We list the explicit forms of the non-zero components of the Lie derivatives C± (along 
the direction of l±) of the 1-form basis: 

{C+df) M = [d M l%) = 0, 

(C.df) M = (d M f_) = ffy-L 0, 0, 0, U -L) , 

(c + du) M = (d M ii) = fa-^, o, o, o, 9*—^) , 

(L„du) M = (ft, C ) = (-*^^, 0, 0, 0, -^^f^) . (156) 

Therefore, 

(At * (3) 1) = [(«*$ + u v ± d v ) e b u 2 {u + f)] rfx 1 A dx 2 A dx 3 . (157) 
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The Hodge dual with respect to h should bqfj 



* (3) 1 = \/det (3) hdx 1 A dx 2 A dx 3 = e b u 2 (u + f)dx l A dx 2 A dx 3 . (158) 

The expansion is given as the *^-dual of the Lie derivatives C± (along the direction 
of l±) of the perpendicular volume form (I158P 



e ± = *® ( £± ,o) ^ = x [(^a. + /« 9u ) e \\ u + f) ] . (159) 
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